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The Pierre Auger observatory has presented evidence that the arrival directions of cosmic rays 
with energies in excess of 6 x 10 7 TeV may be correlated with nearby active galactic nuclei (AGN). 
In this context we revisit a suggestion based on gamma ray observations that nearby Fanaroff- 
Riley I galaxies such as Cen A and M87 are the sources of the local cosmic rays. We compute 
the accompanying neutrino flux and find a flux within reach of second-generation kilometer-scale 
neutrino telescopes. 



I. GAMMA-RAY OBSERVATIONS 

Recently the Pierre Auger observatory has presented 
evidence that the arrival directions of cosmic rays with 
energies in excess of 6 x 10 7 TeV may be correlated with 
nearby active galactic nuclei (AGN) [l|] . The skymap sug- 
gest that events cluster in the direction of the Fanaroff- 
Riley I (FRI) radio galaxy Centaurus A. This source, 
along with the other nearby FRI M87, has previously 
been singled out as a potential cosmic ray accelerator on 
the basis of gamma-ray data [2J . We revisit those argu- 
ments here as well as asking the question whether or not 
kilometer-scale neutrino observatories have the capability 
to unambiguously detect these sources. 

It has been speculated for some time that supermas- 
sive black holes at the centers of active galaxies power 
particle flows that create the opportunity for particle ac- 
celeration to super- EeV energy Q . Both the accretion 
flow near the black hole nucleus and the extended jets 
are potential sites for shocks that allow for efficient par- 
ticle acceleration. Blazars, a subclass of FRI whose jets 
are oriented along our line of sight, are abundant sources 
of TeV energy photons. Exceptionally, the nearby (16 
Mpc) FRI M87 was observed in the late 1990's by the 
HEGRA stereoscopic system of five imaging atmospheric 
Cherenkov telescopes despite the fact that its jet is an- 
gled to more than 30° from our line of sight [4] . The 
integrated flux observed is [f| 
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Recently, the H.E.S.S. array obtained a limit of lO - " 1 m 
the same units on Cen A [6j , a similar source in the south- 
ern hemisphere at a distance of only 3.4 Mpc. Evidence 
for TeV emission from Cen A was obtained in the early 
1970's with the Narrabri optical intensity interferometer 
of the University of Sydney 7] , with an average flux of 



F 1 (E 1 > 0.3 TeV) = 4.4 ± 1.0 x 10 
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The variable flux was collected in two periods of height- 
ened activity lasting roughly one year, pointing at a re- 
gion of coherent emission of size of order 0.3 pc. This 
is consistent with the idea that the high energy emission 
is from an isotropic region near the central black hole 
of mass Mbh = 2 x 10 s solar masses, about two orders 



of magnitude more massive than the one at the center 
of our galaxy ■ The fact that the detector beam did 
not include the radio lobes at the ends of the jets further 
supports the idea of a central engine at the base of the 
jets. The TeV fluxes of Cen A and M87 imply similar 
source luminosities of roughly 7 x 10 40 ergs _1 , assuming 
an E~ 2 gamma-ray spectrum. This suggest that they 
are generic AGN, a fact we will exploit to construct the 
diffuse neutrino flux from all FRI. 



II. MODELING THE SOURCES 

We will assume that the photon flux is of pionic ori- 
gin, i.e. that the high energy protons hinted at by the 
Auger data produce neutral pions on the gas surround- 
ing the central black hole that decay to produce the de- 
tected high-energy photons. This is a plausible model 
given that the proton interaction length on a particle 
density of n ~ 10 6 cm -3 is 4.6 pc, not much larger than 
the source of size R ~ 0.3 pc. We conservatively as- 
sume a density similar to the one near the center of 
our own galaxy. The probability that the accelerated 
proton interacts is therefore expected to be of order 
P ~ R/Xmt ~ Rnap ro t ~ 10%. The energy produced in 
a one year burst is 

-^burst — -PGAfe/f M in faii/ R. 

This leads to M; n f a n ~ 5 x 10 ~ 5 solar masses. This is two 
orders of magnitude larger than the estimated infall at 
the center of our own galaxy [8j and is reasonable because 
of the burst nature and the larger mass of the Cen A black 
hole. We tentatively conclude that the central accelerator 
in FRI such as Cen A and M87 feeds a beam dump and 
pursue the implication of the Auger and TeV gamma 
ray observations for the emission of neutrinos, inevitably 
produced by the charged pions produced in association 
with the neutral ones . Given that the black hole in 
Cen A is more massive than the one at the center of our 
own galaxy by two orders of magnitude, our assumption 
for the target density should therefore be considered a 
lower limit and values of P larger by up to two orders of 
magnitude should be acceptable. 

We will first estimate the neutrino flux from individ- 
ual nearby sources and subsequently derive a diffuse flux 
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assuming that the nearby sources are representative of 
the extragalactic source population. We finally comment 
on the suggestion that the pions are produced in photo- 
hadronic (jyy) rather than hadronic (pp) interactions [lol ]. 

The fluxes of neutrinos and gamma rays of pionic origin 
in astrophysical sources are related to the initial acceler- 
ated proton spectrum via their final-state multiplicities 
and average fractional energies Xj relative to the proton. 
Hadronic collisions leads to equal numbers of each type 
of pion: p + p — > 7r° + ir + + 7r~ . The neutral pion decays 
into two gamma-rays ir° — > 7 + 7 and the charged pions 
into leptons and neutrinos tt^ — > e + 2/^ + + 24. The 
effect of oscillations on the neutrino flux is to equalize the 
number of neutrinos of each flavor arriving at Earth, giv- 
ing one neutrino of each flavor per charged pion. Hence 
we have two-thirds of a charged pion per interacting pro- 
ton and one neutrino (of each flavor) per charged pion; 
and one-third of a neutral pion per interacting proton 
and two photons per neutral pion. Therefore, 
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where R p = 1 — e~ T is the fraction of protons that inter- 
act inside the source where r is the average number of 
interactions per proton, analogous to optical depth. We 
derive the fractional energies from the measured pion in- 
elasticity of approximately 0.2 and from the fact that 
the final-state particles each take an equal fraction of 
the pion energy. This gives us x v ~ 0.25 x n ~ 0.05 and 
x~f ~ 0.5 x-k ~ 0.1. 

To calculate the initial proton flux we must use 
the observed Auger cosmic-ray flux, which is what re- 
mains of the protons after interacting inside the source: 
dN cr /dE = e~ T dN p /dE. We calculate the flux as- 
suming a power-law spectrum of the form dN/dE = 
Nq(E I Eq)" , where we can derive the normalization iVo 
from the observational parameters of the Auger exper- 
iment by calculating the number of events observed [l[ • 
For a — —2.0, we have N cvcnts = field of view x time x 
efficiency x iVo/-Bthresh- This gives us dN cr /dE = 
Advents x 10- 13 (E/TeV)- 2 TcV- 1 cm" 2 s" 1 . Since the 
gamma-ray flux must be below the limit set by HESS, 
we can constrain the unknown source parameters N and 
r: 



introduced. Plugging these constants into the expression 
for the neutrino flux per flavor, we find 
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With a location at the edge of Auger's sky coverage, 
a similar calculation for M87 is at present not possi- 
ble. Clearly, we anticipate the observation with increased 
statistics of events correlated with M87. 

Repeating this calculation for power-law indices be- 
tween -2.0 and -3.0, we obtain between 0.8 and 0.02 
events/year for a generic neutrino detector of effective 
muon area ~ 1 km 2 . 



III. DIFFUSE FLUX 

Having determined the neutrino flux from a UHE cos- 
mic ray source, we can determine the total diffuse flux 
from all such sources within the horizon. We assume here 
that the common luminosities of Cen A and M87 are rep- 
resentative of FR1 radio galaxies. Given an FRI de nsity 
of n ~ 8 x 10 4 Gpc -3 within a horizon of R ~ 3 Gpc pl| . 
the total diffuse flux from all 4-7rsr of the sky is simply 
the sum of the luminosities of the sources weighted by 
their distances: 
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where we perform the sum by assuming the galaxies are 
uniformly distributed inside the sphere of the horizon. 
For a Centaurus neutrino flux that goes as E~ 2 , this 
evaluates to 
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approximately one-tenth of the Waxman-Bahcall 
flux [l2l|. Varying the spectral indices as before, we 
obtain an event rate per km 2 year from the northern sky 
of between 19 and 0.5 neutrinos (Figure [T|). Considering 
sources out to 3 Gpc, a redshift of order 0.5, is probably 
conservative and extending the sources beyond z ~ 1 
may increase the flux by a factor 3. 



IV. PHOTOHADRONIC PROCESSES 
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This gives us the constraint N(e T — 1) < 150. The 
number of events correlating to Cen A depends on the 
unknown intergalactic magnetic fields. Two correlate 
within 3 degrees but as many as ten could be associated 
with the source. In this case r should be C(a few) which 
is consistent with the physics of the problem previously 



Often analyses of AGNs assume that any production of 
secondary particles that serve as a possible signal comes 
from photohadronic interactions of accelerated charged 
particles with ambient radiation fields. This means that 
it is possible to model the emission, in both gamma-rays 
and neutrinos, of an optically thick ultra-high-energy cos- 
mic ray source Q . 

Deriving the neutrino and gamma-ray spectra from the 
proton spectrum proceeds as in the p-p case, with the dif- 
ference that we must take the different branching ratios 
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FIG. 1: Total diffuse neutrino flux and event rate from FRI's 



of the delta channel p + 7 — > A + — > + 7r°(7r + ) into 
account. Due to isospin conservation, the 7r° final state 
is twice as likely as the n + (the opposite of the p-p case), 
resulting in a relative factor of 4 [13[ . 

One general feature of such a process is that we ex- 
pect breaks in the cosmic-ray spectrum corresponding 
to the onset of energy loss processes, in other words, 
when the cosmic ray primary energy is large enough to 
produce pions on the photon field. Following Ref. [Ic| |. 
we can then model the gamma-ray and neutrino spectra 
with a break at the corresponding secondary energy. If 
we vary the spectral index of the protons, assuming the 
index hardens by one for energies below the break, we 
can obtain the neutrino spectrum by fixing the cosmic- 
ray spectrum above the break at the Auger measurement 
and fixing the gamma-ray spectrum below the break to 
the Narrabri measurement [3] . If we take the break en- 
ergy as variable within a decade of 10 5 TeV for protons 
(~ 10 4 TeV for gamma-rays) and allow that the actual 
number of cosmic-ray events from Cen A is likely different 
from the Auger experiment's estimate due to the strength 
of the intergalactic magnetic field being unknown, we can 
ask which combinations of spectral index, break energy 
and number of cosmic ray events leads to a continuous 
spectrum of neutrinos and gamma rays. The final an- 
swer is relatively insensitive to the number of cosmic ray 
events and most sensitive to the spectral index, allow- 



ing us to get a spectrum by picking a spectral index and 
then tuning with the break energy and cosmic ray events. 
The neutrino event rate per neutrino flavor per km 2 year 
increases with harder spectra and lower break energies, 
giving ~ 5 events/year for a proton spectrum with a a 
break at 4 x 10 3 TeV and a spectral index of -2.7 above 
the break, and ~ 0.03 events/year for a proton spectrum 
with a break at 5 x 10 5 TeV and a spectral index of - 
3.0 above the break. We can take these event rates as 
extreme since the break energies required for continuity 
are outside the range deemed reasonable. Proton spec- 
tral indices harder than -2.7 above the break do not give 
continuous spectra for break energies above 10 3 TeV and 
so were disregarded. 

This analysis may be deficient in several regards. On 
the technical side, it does not account for the cosmic 
rays being the remnant of the proton spectrum that did 
not interact in the source since it does not attempt to 
model the optical depth of the source. Moreover, if the 
spectral break is due to the onset of pionization, there 
should be no gamma-rays with energies below x-yEb rea k, 
as the source is transparent to the protons at low ener- 
gies. Thus, in this model the observed TeV gamma ray 
fluxes from M87 and Cen A cannot be explained. More 
sensitive observations of Cen A by H.E.S.S. are awaited 
and confirmation of a TeV flux at the levels assumed in 
this paper should disfavor the photoproduction model. 

In summary, we conclude that FRI modeled on the 
fluxes measured by the Auger and TeV gamma ray 
telescopes predict neutrinos at the level of 10% of the 
Waxman-Bahcall flux, possibly higher by the factor re- 
lating to cosmic evolution mentioned earlier in this paper. 
There should be no doubt that the predicted flux should 
be within reach of IceCube [3| and a future Mediter- 
ranean kilometer-scale neutrino telescope, as it is at the 
level of IceCube's 90% confidence level sensitivity in a 
single year . It follows that detection at the 5<r level 
should occur within ~ 5 years. 
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